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Abstract 

It was recently proposed that the Euclidian norm of the tj vector of the 
coupled cluster wave function (normalized by the number of electrons included in 
the correlation procedure) could be used to determine whether a single-reference- 
based electron correlation procedure is appropriate. This diagnostic, Tj, is defined 
for use with self-consistent-field molecular orbitals and is invariant to the same 
orbital rotations as the coupled cluster energy. Tj is investigated for several different 
chemical systems which exhibit a range of multireference behavior, and is shown 
to be an excellent measure of the importance of non-dynamical electron correlation 
and is far superior to C 0 from a singles and doubles configuration interaction wave 
function. It is further suggested that when the aim is to recover a large fraction of 
the dynamical electron correlation energy, a large T) (be., > 0.02) probably indicates 
the need for a multireference electron correlation procedure. 
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Introduction 


It was recently proposed 3 that the Euclidian norm of the vector of tj ampli- 
tudes in the closed-shell coupled cluster singles and doubles wave function could be 
used as a diagnostic for the a priori prediction of the reliability of results obtained 
from a single-reference-based electron correlation procedure. The tj amplitudes in 
coupled cluster theory are closely related to the coefficients of singly excited config- 
urations in configuration interaction theory. It is well documented 2 that the singly 
excited configurations in an electron correlation procedure allow molecular orbital 
relaxation to occur. For many years quantum chemists have used Co, the reference 
configuration coefficient in a configuration interaction wave function, as a diagnos- 
tic. As is widely recognized, however, if Co is taken from a self-consistent-field 
(SCF) singles and doubles configuration interaction (CISD) wave function, then 
it is of limited utility since the molecular orbitals are strongly biased towards the 
SCF reference function. Thus it is not uncommon for a known multireference sys- 
tem to yield an SCF-CISD Co which is 0.95 or larger (i.e., the SCF determinant 
comprises 90% of the wave function). A reliable diagnostic which is more sensitive 
to the importance of non-dynamical electron correlation would therefore be of great 
utility. 

Laidig, Purvis and Bartlett 3,4 have investigated the use of localized molecular 
orbitals in coupled cluster methods, specifically the doubles and the singles and 
doubles coupled cluster methods (CCD and CCSD, respectively). We note that 
the particular localization technique investigated by Laidig et aL does not leave 
the SCF energy unaffected 3,4 . They found that the use of localized molecular 
orbitals greatly improved the CCD results, but that the CCSD energies were little 
affected by the different reference molecular orbitals. The inclusion of e Tl in the 
CCSD wave function thus accounts for the important orbital relaxation effects which 
were incorporated by localizing the molecular orbitals. In addition, Scuseria and 
Schaefer^ have investigated the use of Brueckner-like molecular orbitals in CCSD 
and CCD calculations and arrive at essentially the same conclusions. 

The purpose of the present study is to further investigate the use of the Euclid- 
ian norm of t } as a diagnostic: applying this test to chemical systems exhibiting a 
range of bonding situations and known multireference and strongly single-reference 
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dominated problems. In tins way the actual value and utility of the Euclidian norm 
of tj as a diagnostic tool will become evident. To begin, it must be emphasized 
that the diagnostic reported here 


was always determined using SCF molecular orbitals. As the results of references 
3 through 5 clearly demonstrate, it is possible to obtain a similar C-CSD energy 
with different molecular orbitals which will give a different tj vector and a different 
Euclidian norm. In fact, the Euclidian norm of tj for the ‘optimized orbitals’ 
of reference 5 should be very' close to zero. Thus, in order to compare Tj from 
different chemical systems, the diagnostic must be uniquely defined for each system. 
The most straightforward. approach is to require that restricted Hartree-Fock SCF 
molecular orbitals are used to determine the CCSD wave function and diagnostic 
for each system, and this is therefore the approach which has been adopted in the 
present study. In addition, we point out that since the CCSD energy is invariant to 
unitary transformations of occupied-occupied or virtual-virtual molecular orbitals 
the Ti diagnostic will also be invariant to these types of orbital rotations. 

The next section contains a brief summary of the theoretical methods used 
together with a description of the method we have devised to judge the 7j diag- 
nostic. The results, including a discussion, are presented m the third section. Our 
conclusions are presented in the final section. 

Methods 

All of the chemical systems included in this study have been investigated 
previously 1,6 ' 7 and these reports include a detailed description of the basis sets 
and geometries. We therefore include only a brief description of the basis sets. 
Table 1 contains the size of the primitive basis, our designation and the reference 
from which the orbital exponents and contraction coefficients may be obtained. 
In forming the designation for each basis two rules have been followed. Firstly, a 
generally contracted atomic natural orbital (ANO) basis set is denoted by square 
brackets, e.g., [4321], where the numbers enumerate the number of contracted s, p, 
d and / functions, respectively. Secondly, a basis set which utilizes a segmented 
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contraction scheme is designated as 7s3p2dlf. for example. In most cases where a 
segmented contraction is used, the contraction has been performed over the core 
atomic orbitals, allowing maximum flexibility in the valence region. For those cases 
where the polarization function orbital exponents are not given in the reference 
the exponents are listed in Table 1. In addition, where more than one level of 
polarization function has been included (e.g., 7s3p2dlf Be) the levels are separated 
by a semicolon. 

Bond lengths are given in atomic units, a 0 . The unique bond length is specified 
for the Be 3 , Mg 3 , Be 4 and Mg 4 clusters. The trimers form an equilateral triangle 
and the tetramers adopt a tetrahedral structure. The pentamer, Be 5 , is defined by 
two bond lengths since it conforms to a trigonal bipvramidal geometry. For this 
system, the first bond length refers to a side of the triangular base while the second 
refers to the distance from an apex atom to one contained in the base. The bond 
lengths and bond angles for FOOF, (N0)2 and FNNF are the TZ2P MP2 structures 
reported in reference 1. 

Since the definition of depends upon the number of electrons correlated 
it is clearly important to consider which electrons should be included in this def- 
inition. It is expected that only the valence electrons should be important for 
non-dynamical electron correlation effects and therefore we have chosen to freeze 
the core-like molecular orbitals in all procedures. It is possible that even if the core 
electrons are included in the correlation procedure, then the definition of T\ should 
include only the number of valence electrons (see note added in proof in reference 1). 
However, as several studies have demonstrated 218,9 the basis set requirements for the 
adequate treatment of core-valence and core-core correlation effects are quite severe. 
Therefore, for our initial investigations of T\ only the valence electrons are consid- 
ered. Additionally, for basis sets which utilize segemented contractions the virtual 
molecular orbitals which are the core-counterparts were deleted from the correlation 
procedure. The CCSD wave functions were determined with a vectorized closed- 
shell CCSD method 10 and the Cl wave functions were evaluated with either the 
Berkeley shape-driven graphical unitary group Cl program 11 or the MOLECULE- 
SWEDEN codes 12,13 . 
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Results and Discussion 


The Ti diagnostic together with Co from CISD and full Cl wave functions are 
presented for several systems in Table 2. Note that only two electrons are correlated 
for the first five systems. Comparing Tj with Co for these systems it is clear that 
there is a good correspondence between and the total weight of the reference in 
the full Cl wave function. Thus for He and H 2 , where C 0 is greater than 0.99, the 
Tj diagnostic is 0.0029 and 0.0050, respectively, whereas for the other three systems 
(Be, Mg and Li 2 ) X 3 is greater than 0.015 and C 0 is less than 0.965. Be and Mg are 
known to exhibit multireference behaviour due to the s — p near degeneracy. Li 2 
possesses a a — a* near degeneracy in addition to the s — p near degeneracy. 

Since the remaining molecules in Table 2 (He 2 , Be 2 , Mg 2 and HF) all contain 
more than two valence electrons it is possible to compare Tj , the Co from CISD, and 
the C 0 from full Cl. As expected, He 2 is strongly dominated by a single reference 
function and so there is not a significant difference between the CISD and full Cl 
Cp. Consistently, T\ is again very small and is actually the same (to the precision 
reported) as for the single He atom. However, for the Be 2 , Mg 2 and HF diatomics 
there is a significant difference between the full Cl and CISD Co* In fact, for 
Be 2 the difference amounts to 4.5% of the full Cl wave function. An important 
point which should be emphasized is that due to the lack of size-extensivity the 
discrepancy between the full Cl and CISD Co is expected to become larger as the 
number of electrons correlated increases. is greater than 0.013 for Be 2 , Mg 2 and 
HF demonstrating that a large degree of orbital relaxation occurs. Thus, the results 
of Table 2 demonstrate two important points: 1) there is a good correspondence 
between X 3 and the full Cl Co when a modest number of electrons are correlated 
and 2) for chemical systems with more than two electrons there may be a large 
difference between the CISD and full Cl Co. 

T\ and the CISD Co for several different chemical systems with a large number 
of valence electrons are collected in Table 3. The T 3 diagnostic and the Co for the 
Be and Mg clusters (at their equilibrium structures) indicates that these systems 
are probably not well described by a single-reference method and that a large degree 
of orbital relaxation is taking place. Binding energies and equilibrium bond lengths 
for the clusters, for example, would be expected to be substantially in error when 
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a single-reference-based treatment is used, and it is doubtful that binding energy 
predictions would be reliable to within even 10 kcal/mol. However, the large Co 
for these systems might tempt many observers to believe that a single-reference- 
based electron correlation procedure is adequate. Conversely, the T\ value is larger 
than 0.02 for each cluster with the exception of Mg 3 . The comparisons made in 
Table 2 together with multireference Cl (MRCI) results 6 suggest that multireference 
techniques are required for Be and Mg clusters, and thus that a Tj value larger 
than 0.02 is a clear indication that other import ant configurations exist and may be 
needed as references in a treatment of dynamical electron correlation. The infinite 
separation results for the Be and Mg clusters also demonstrate the inadequacy of 
using C 0 from a CISD wave function since the C 0 suggests that as the number 
of atoms increases the ‘super-molecule’ is more difficult to describe whereas the 
size-extensive CCSD method correctly shows that these systems are all equivalently 
described (in fact, since only valence electrons are correlated the CCSD results 
correspond to a full Cl). 

The FOOF, (NO) 2 and FNNF molecules are included in Table 3 since these 
were the systems investigated in the study 1 which first suggested the use of Tj . 
These systems are very difficult to describe — the geometry of FOOF is not even 
qualitatively correct at the CISD level, for example. The C 0 values for FOOF and 
(NO) 2 are very similar although Tj indicates that non-dynamical electron corre- 
lation is much more important for FOOF. The results of several single-reference 
methods for these two systems 1 provide additional evidence that the electron cor- 
relation of FOOF is indeed even more difficult to describe than that for (NO) 2 . 
The C 0 for the isomers of FNNF suggests that these systems are more strongly 
dominated by a single reference than either FOOF or (NO) 2 and that they are all 
nearly equally well described by a single-reference-based method. However, while 
the Tj diagnostic does suggest that non-dynamical electron correlation is less im- 
portant in the cis and trans isomers, it also indicates that the transition state is 
strongly affected by non-dynamical electron correlation and thus, single-reference- 
based methods will not work as well for TS-FNNF as they do for cis and trans 
FNNF. Again, the latter conclusions are consistent with the results of reference 1. 

The last four molecules of Table 3 are all known to be strongly dominated 
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by a single determinant reference function and both C 0 and Tj are consistent with 
this observation. However, the fine details of relating Co and Tj exhibit small 
inconsistencies. For example, for the first-row closed-shell hydrides it is generally 
accepted that the reliability of a single-reference-based electron correlation method 
decreases in the order CH4 > H 2 0 > HF. The T\ diagnostic is consistent with this 
empirical observation whereas the C 0 from a CISD wave function exhibits exactly 
the opposite trend. 

Perhaps the molecule which best exhibits the superiority of T\ over Co as a 
diagnostic is the CuH diatomic. It has been shown 14 that the bonding of CuH is 
complicated because of the importance of both the d? s 1 and d 10 ^ 1 atomic occu- 
pations of Cu. Thus there are several important configurations which differ from 
the closed shell single determinant reference by a single excitation. The C 0 for 
this diatomic is 0.96, w’hich is very similar to that obtained for CH 4 . However, 
the Ti value, 0.046, is the largest found in this study. Thus, the important non- 
dynamical electron correlation effects present in the bonding of CuH are completely 
missed by the single-reference CISD method whereas the T x diagnostic correctly 
indicates the importance of these effects. 

Table 4 contains T\ and Co (from CISD) for Be 3 and H 2 0 using several basis 
sets in order to determine the one-particle basis set effect. In order for this diagnostic 
to be generally useful it should exhibit a certain degree of invariance with respect 
to the choice of a one-particle basis set. This statement assumes, of course, that 
the smallest basis sets at least contain proper correlating functions. On the other 
hand, it is well known that the one-particle and n-particle basis sets are inherently 
coupled though this coupling is usually small. As the one-particle basis set limit 
is approached, it may be expected that T\ will stabilize. This is expected despite 
the fact that the n-particle basis increases substantially with one-particle basis set 
augmentations. 

The results of Table 4 confirm the above discussion and demonstrate that 
T] converges to a value near 0.0340 for Be 3 and near 0.0075 for H 2 0. The fact 
that Ti decreases with improvements in the one-particle basis set provides further 
support for the above discussion. In other w T ords, a larger degree of orbital relaxation 
is required for the smaller one-particle basis sets (giving a larger T ) in order to 


/ 


compensate for the lack of flexibility. Thus the value to which X 3 converges should 
give an indication of the inherent importance of non-dvnamical electron correlation 
for the chemical system under investigation. Moreover, the rate of convergence with 
respect to basis set improvement should give a measure of the interaction between 
the one- and n-particle basis sets. 

Conclusions 

The Xi diagnostic has been shown to be a reliable measure of the impor- 
tance of non- dynamical electron correlation and to be far superior to the use of Co 
from a CISD wave function as an indicator as to whether it is appropriate to use 
a single-reference-based electron correlation procedure. No doubt a similar type of 
diagnostic could be defined for the CISD wave function by separating the C 3 coeffi- 
cients (coefficients from the singly excited configurations). However, there are two 
problems with this procedure. Because CISD is, in general, not size extensive then 
the diagnostic would not have the desired property of giving the same result for two 
non-interacting He atoms as it would for a single He atom. Also, the results for 
CuH presented in this study clearly indicate that the CISD procedure is incapable 
of overcoming the bias of using SCF molecular orbitals and thus any diagnostic sim- 
ilar to X 3 but based on an SCF-CISD wave function would almost certainly suffer 
from this bias. 

Several studies 14,22 have pointed out that the coupled pair functional 23 (CPF), 
modified CPF 24 (MCPF) and averaged CPF 25 (ACPF) methods are very good at 
identifying specific configurations which are important and hence should be used as 
references in a multireference electron correlation procedure. A diagnostic similar 
to A, could also be constructed for these methods and it is likely that it would give 
similar results to X for those situations where the CCSD and the various CPF- 
tvpe methods gave similar results. Clearly for the situation where the methods 
give verv different results (such as FOOF 1 ) the diagnostics would be expected 
to yield different results also. In addition these observations also indicate that 
specific important configurations may be identified by analysis of the CCSD 1 3 and 
i 2 amplitudes. 

Finally, the results of this study indicate that if X, is greater than 0.02 then 
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single-reference-based electron correlation methods are probably unreliable and will 
certainly not yield highly accurate results. 

Acknowledgements . 

PRT was supported by NASA grant NCC2-371 and TJL was partially sup- 
ported by NASA grant NCC2-552. We would like to thank the NAS Facility for 
early access to the CRAY Y-MP. 


9 


References 


1. T. J. Lee, J. E. Rice, G. E. Scuseria and H. F. Schaefer, Theor. Chim. Acta. 
75, 81 (1989). 

2. I. Shavitt in “Methods of Electronic Structure Theory,” Vol. 3, Editor H. F. 
Schaefer (Plenum. New \'ork, 1977). 

3. W. D. Laidig, G. D. Purvis and R. J. Bartlett, International J. Ouan. Chem. 
Symp. 6, 561 (1982). 

4. W. D. Laidig, G. D. Purvis and R. J. Bartlett, Chem. Phys. Lett. 97, 209 
(1983). 

5. G. E. Scuseria and H. F. Schaefer, Chem. Phys. Lett. 142, 354 (1987). 

6. T. J. Lee, A. P. Rendell and P. R. Taylor, to be published. 

7. C. W. Bauschlicher, S. R. Langhoff, H. Partridge and L. A. Barnes, submitted 
to J. Chem. Phys. 

8. H. Partridge, C. W. Bauschlicher, S. P. Walch and B. Liu, J. Chem. Phys. 
79, 1866 (1983). 

9. C. W. Bauschlicher, S. R. Langhoff and P. R. Taylor, J. Chem. Phys. 88, 
2540 (1988). 

10. T. J. Lee and J. E. Rice, Chem. Phys. Lett. 150, 406 (1988). 

11. P. Saxe, D. J. Fox. H. F. Schaefer and N. C. Handy, J. Chem. Phys. 77, 5584 
(1982). 

12. J. Almlof, MOLECULE, a vectorized Gaussian integral program. 

13. SWEDEN is a vectorized SCF-MCSCF-direct CI- conventional CI-CPF- 
MCPF program written by P. E. M. Siegbahn, C. W. Bauschlicher, B. Roos, 
P. R. Taylor, A. Heiberg, J. Almof, S. R. Langhoff and D. P. Chong. 

14. D. P. Chong, S. R. Langhoff, C. W. Bauschlicher, S. P. Walch and H. Par- 
tridge. J. Chem. Phys. 85, 2850 (1986). 

15. S. Huzinaga, J . Chem. Phys. 42, 1293 (1965). 

16. T. Dunning, J. Chem. Phys. 53, 2823 (1970). 

17. F. B. van Duijneveldt., IBM Res. Dept. RJ 1971, 945. 

18. J. Almlof and P. R. Taylor, J. Chem. Phys. 86, 4070 (1987). 

19. T. Dunning, J. Chem. Phys. 55, 716 (1971). 


10 



20. A. D. McLean and G. S. Chandler. J. Chem. Phys. 72, 5639 (1980). 

21. C. W. Bauschlicher, S. R. Langhoff, P. R. Taylor, N. C. Handy and P. J 
Knowles, J. Chem. Phys. 85, 1469 (1986). 

22. C. W. Bauschlicher, S. R. Langhoff, T. J. Lee and P. R. Taylor, J. Chem 
Phys., in press. 

23. R. Aldrichs, P. Scharf and C. Ehrhardt, J. Chem. Phys. 82, 890 (1985). 

24. D. P. Chong and S. R. Langhoff, J. Chem. Phys. 84, 5606 (1986). 

25. R. J. Gdanitz and R. Aldrichs, Chem. Phys. Lett. 143, 413 (1988). 


11 


Table 1 


Basis set designations and definitions used in this study. 


Atom 

Primitive Basis 

Designation 

Reference 

Polarization 

Exponents 

H 

4s 

DZP 

15.16 

0.75 

H 

8s2p 

6s2p 

17 

1.0,0.33 

H 

8s6p 

[32] 

18 

- 

H 

8s6p4d 

[321] 

18 

- 

H 

8s6p4d 

[432] 

18 

- 

He 

8s2p 

6s2p 

17 

1.0,0.33 

Li 

9s4p 

4s3p 

17 

- 

Be 

12s5p2d 

7s3p2d 

17 

0.3, 0.1 

Be 

12s5p2dlf 

7s3p2dlf 

17 

0.3,0.1;0.26 

Be 

12s7p4d2f 

[5321] 

18 

- 

Be 

12s7p4d2f 

[6432] 

18 

- 

C 

10s6p2d 

TZ2P° 

15.19 

1.5,0.35 

X 

10s6p2d 

TZ2P° 

15.19 

1.5,0.35 

0 

10s6p2d 

TZ2P° 

15,19 

1.5,0.35 

0 

13s8p6d 

[432] 

18 

- 

0 

13s8p6d4f 

[4321] 

18 

- 

0 

13s8p6d4f2g 

[54321] 

18 

- 

F 

9s5pld 

DZP 

15,16 

1.6 

F 

10s6p2d 

TZ2P° 

15,19 

1.5.0.35 

Ne 

10s6p2d 

TZ2P 6 

17 

4. 5, 1.3 

Mg 

12s9p2d 

6s5p2d 

20 

0.3, 0.1 

Cu 

14sll P 6d3f 

8s6p4dlf 

7 



a. The 5s3p contraction of reference 19 was used. 

b. A 5s3p contraction, similar to those given in reference 19, was constructed. 
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Table 2 

The Xi diagnostic together with the Co obtained 
from a full Cl and a CISD wave function 0 . 


Molecule 

Basis Set 

T 

T, 

c b 

cs 

He 

6s2p 

- 

0.0029 

0.9960 

0.9960 

Ho 

6s2p 

1.361 

0.0050 

0.9912 

0.9912 

Be 

7s3p2d 

- 

0.0210 

0.9523 

0.9523 

Mg 

6s5p2d 

- 

0.0159 

0.9640 

0.9640 

Ll o 

4s3p 

5.11 

0.0165 

0.9510 

0.9510 

He 2 

6s2p 

5.61 

0.0029 

0.9920 

0.9921 

Be 2 

7s3p2d 

4.75 

0.0282 

0.8901 

0.9150 

Mgo 

6s5p2d 

7.35 

0.0138 

0.9268 

0.9401 

EF d 

DZP 

2.5995 

0.0187 

0.9583 

0.9680 


a. All correlated wave functions are based upon SCF molecular orbitals. Only 
valence electrons have been included in the correlation procedure. Bond lengths are 
in atomic units, a- 0 . 

b. FCI. 

c. CISD. 

d. FCI and CISD results from reference 21. 
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Table 3 

The T\ diagnostic together with C 0 obtained from a C1SD wave function 0 . 


Molecule 

Basis Set 

Geometry 

T, 

Co 

CuH 

b 

2.850 

0.0461 

0.9621 

Be 3 

7s3p2d 

4.273 

0.0360 

0.9133 

Be 4 

7s3p2d 

3.915 

0.0318 

0.9189 

Be 5 

7s3p2d 

3.831,3.929 

0.0290 

0.9094 

Be 3 

7s3p2d 

OG 

0.0210 

0.9067 

Be 4 

7s3p2d 

oo 

0.0210 

0.8933 

Be 3 

7s3p2d 

oo 

0.0210 

0.8828 

Mg 3 

6s5p2d 

7.522 

0.0127 

0.9235 

Mg 4 

6s5p2d 

6.102 

0.0204 

0.9102 

Mgs 

6s5p2d 

OG 

0.0159 

0.9240 

Mg4 

6s5p2d 

OC 

0.0159 

0.9111 

FOOF 

TZ2P 

see text 

0.0313 

0.9189 

(NO ) 2 

TZ2P 

see text 

0.0203 

0.9177 

cis-FNNF 

TZ2P 

see text 

0.0187 

0.9303 

trans-FNNF 

TZ2P 

see text 

0.0166 

0.9308 

TS-FNNF 0 

TZ2P 

see text 

0.0277 

0.9283 

HF 

TZ2P 

1.734 

0.0104 

0.9775 

H 2 0 

TZ2P 

1.809,104.8° 

0.0096 

0.9720 

ch 4 

TZ2P 

2.052 

0.0073 

0.9672 

Ne 

TZ2P 

- 

0.0065 

0.9850 


a. All correlated wave functions are based upon SCF molecular orbitals. Only 
valence electrons have been included in the correlation procedure. Bond lengths are 
in atomic units, ao- 

b. The Cu basis is as described in Table 1 and the H basis is the [32] ANO basis 
set. 

c. Transition state to cis-trans isomerization. 
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Table 4 

The T] diagnostic for Be 3 and H 2 O using several different 
basis sets. The Cq value is obtained from a CISD wave function. 


Molecule 

Basis Set 

T, 

Co 

Be 3 

7s3p2d 

0.0360 

0.9133 

Be 3 

7s3p2dlf 

0.0339 

0.9149 

Be 3 

[421] 

0.0386 

0.9107 

Be 3 

[5321] 

0.0341 

0.9148 

Be 3 

[6432] 

0.0341 

0.9157 

H 2 0 

TZ2P 

0.0096 

0.9720 

h 2 o 

[432/32] 

0.0076 

0.9721 

h 2 o 

[4321/321] 

0.0071 

0.9714 

h 2 o 

[54321/432] 

0.0078 

0.9713 


a. The geometries are the same as those listed in Table 2. Only valence electrons 
have been correlated. 
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